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Abstract We obtain improved constraints on the cou¬ 
pling constants of axion-like particles to nucleons from 
a recently performed Casimir-less experiment. For this 
purpose, the differential force between a Au-coated sphe¬ 
re and either Au or Si sectors of a rotating disc, arising 
due to two-axion exchange, is calculated. Over a wide 
region of axion masses from 1.7 x 10“^ eV to 0.9 eV 
the obtained constraints are stronger up to a factor of 
60 than the previously known ones following from the 
Cavendish-type experiment and measurements of the 
effective Casimir pressure. 


1 Introduction 

It is common knowledge that the proper QCD axions 
are pseudoscalar particles which appear as a conse¬ 
quence of breaking the Peccei-Quinn symmetry [T] pro¬ 
posed to resolve the problem of strong CP violation 
in QCD. After the prediction of axions m. a lot of 
experimental and theoretical work has been done on 
their search and investigation of their role in elementary 
particle physics, astrophysics and cosmology [lll21ini[71ISl 
At the moment, the originally introduced 
QCD axions, which are pseudo-Nambu-Goldstone bosons, 
are constrained to a very narrow band in parameter 
space and many types of so-called axion-like parti¬ 
cles are proposed in different models (see, for instance, 
the hadronic axions [Hill] and the GUT axions [13 

US]). 

Axion-like particles interact with photons, electrons 
and nucleons. Many searches of these particles are based 
on the use of helioscopes and haloscopes [T7] . The helio¬ 
scopes are created for registration of axion-like particles 
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generated in the sun |181I191[^I21) . The haloscopes ex¬ 
ploit an idea that axion-like particles are possible con¬ 
stituents of dark matter wm and fill all the space 
around us. Then, their coupling to photons can be de¬ 
tected using a cryogenic microwave cavity in strong 
magnetic field uimm- Many constraints on the pa¬ 
rameters of axion-like particles were obtained also from 
different astrophysical processes (see, for instance, pSl 
Thus, from the neutrino data 
of supernova SN 1987A the coupling constant of the 
hadronic axions to nucleons was shown to be less than 
10“^° or larger than 10“^ with a narrow allowed region 
in the vicinity of 10“® [26]. Stellar cooling by the emis¬ 
sion of hadronic axions leads to a conclusion that for 
hadronic axions this interaction constant is less than 
3 X 10“^° |271I32|. It is noted however, that the 
emission rate suffers from significant uncertainties re¬ 
lated to dense nuclear matter effects. New strong limit 
on an axion mass and, thus, on an axion-to-nucleon in¬ 
teraction constant (which are connected for hadronic 
axions) was obtained from direct Chandra observations 
of the surface temperature of isolated neutron star in 
Cassiopeia A and its cooling scenario [55] - 

The model-independent laboratory constraints on 
the coupling constants of axion-like particles with nu¬ 
cleons were obtained from neutron physics pi[55] .E6t- 
vos- and Cavendish-type experiments and 

from measurements of the Casimir and Casimir-Polder 
force These constraints cover a wide range 

of masses of axion-like particles from I0“^° eV to 20 eV. 
As was shown in IriiriTTmrn!] (see also [33 for a re¬ 
view) , measurements of the Casimir interaction lead to 
stronger constraints on the coupling constants of axions 
to nucleons than those obtained from the Cavendish- 
type experiments. This corresponds to separation dis- 
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tances between the test bodies where the Casimir in¬ 
teraction becomes stronger than the gravitational one. 

In this paper, we obtain improved constraints on the 
coupling constants of axion-like particles to a proton 
and a neutron following from the recently performed 
Casimir-less experiment [33] . This is the differential force 
measurement between a Au-coated sphere and either 
a Au sector or a Si sector of the structured disc de¬ 
posited on a Si substrate and covered by the overlayers 
of Cr and Au. In such a manner, the contribution of 
the Casimir force to the differential signal is subtracted, 
and the measurement result is determined solely by a 
difference in the forces due to exchange of some hypo¬ 
thetical particles. By achieving the unprecedented force 
sensitivity of approximately 10“^® N, an improvement 
of constraints on the Yukawa-type corrections to New¬ 
tonian gravitation by a factor of 10^ was achieved m- 
The corrections of Yukawa type arise due to exchange 
of one scalar boson between two atoms of the labora¬ 
tory test bodies |3S] or from compact extra dimensions 
with low-energy compactification scale [46] . Here, we 
use the same experimental results to improve the pre¬ 
viously known laboratory constraints on the coupling 
constants of axion-like particles to nucleons. Taking into 
account that the test bodies are unpolarized m and 
the axion-like particles are pseudo-scalar, the additional 
axionic interaction arises due to two-axion exchange be¬ 
tween nucleons of the test bodies. Here, we strengthen 
the axion-to-nucleon coupling constants up to a fac¬ 
tor of 60 within the wide region of axion masses from 
1.7 X 10“^ to 0.9 eV. All equations are written in the 
system of units with h = c = 1. 


2 Differential force between a sphere and a 
structured disc due to two-axion exchange 

In the experiment im, a Au-coated sapphire sphere of 
R = 149.3 fim radius interacts in vacuum with either a 
Au sector or a Si sector of the structured rotating disc 
of thickness D = 2.1 /rm, and the difference in these in¬ 
teraction forces is an immediately measured quantity. 
The structured disc was deposited on a Si substrate 
and covered by the overlayers of Cr and Au of thick¬ 
nesses dcr = 10 nm and (Iau = 150 nm, respectively. 
Note that in [33] the concentric alternating strips of Au 
and Si have been used rather than a sectoral structure. 
This, however, does not influence our calculation of the 
differential force. An important point is that the thick 
overlayer of Au results in equal Casimir forces when the 
sphere bottom is above a Au sector or a Si sector. Thus, 
the Casimir force does not contribute to the measured 
differential force. The latter is determined by possible 


hypothetical interactions, such as the Yukawa-type cor¬ 
rection to Newton’s gravitational law m or the two- 
axion exchange between nucleons of the sphere and the 
structured disc under consideration here (if both these 
attractive interactions exist in nature and contribute 
to the measured differential force, the constraints im¬ 
posed on each of them by the measurement data would 
be even stronger than those obtained in [33] and in this 
paper). 


In this section, we consider the homogeneous Au 
sphere interacting due to two-axion exchange between 
nucleons with the structured Au/Si disc. We assume the 
pseudoscalar character of the axion-nucleon interaction, 
which is applicable to wide classes of axion-like parti¬ 
cles, specifically, to all GUT axions [32] with no con¬ 
nection between their mass and their interaction con¬ 
stant. Note that the account of scalar coupling of axions 
to fermions m or interaction of axions with electrons 
could only slightly increase the magnitude of a differ¬ 
ential axionic force and, thus, only slightly strengthen 
the obtained constraints (see [551H5] for the constraints 
on scalar interaction of axions with nucleons). 


We assume that the coordinate plane (x, y) coin¬ 
cides with the upper plane of the disc and the z axis is 
perpendicular to it. The origin of the coordinate system 
is chosen below the bottom point of the sphere nearest 
to the disc. Without loss of accuracy one can neglect by 
the finite size effects and consider the disc of infinitely 
large area The separation distance between the 

sphere and the disk is a, so that the sphere center is 
at z = a + R. The effective potential due to two-axion 
exchange between two nucleons situated at the points 
ri of the sphere and r 2 of the disc is given by iMiisni 

ED 


Vki{\ri - r-2|) 


glkdli mg Ki{2ma\ri - r 2 |) 
32Tr^m‘^ {ri — r2)^ 


( 1 ) 


Here, gak and gai are the coupling constants of inter¬ 
action between an axion-like particle of mass TOq and 
a proton (fc, I = p) or a neutron (fc, I = n), the mean 
mass of a nucleon is m = {mn + mp)/2, and A'i(z) is the 
modified Bessel function of the second kind. Equation 
ID) is applicable under the condition |ri — r 2 ] ^ 1/m 
which is satisfied with large safety margin because in 
the experiment a > 200 nm [33] • 


The additional force due to two-axion exchange, act¬ 
ing between a homogeneous sphere (s) and a homoge¬ 
neous disc (d), was found in [40] by the summation of 
microscopic forces determined by the potential ID: 

TTTT? /•2R+a 

i^add(a) = / dzi[i?2 - (zi - i? - af] 

xI{ma,D,zi), ( 2 ) 
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where 


d 


I{ma,D,Zi) = — 
uz 


dz-2 


pdp 


1 J-D 


(3) 


Ki{2may/p'^ + (zi - Z2)^) 


+ (zi - Z 2)2 

Here, the coefficient Cd,s for the disc and sphere mate¬ 
rials is defined as 


Cd,s — Pd,s 


9ap ^d,t 


dan 


(4) 


47r pd,s 47r pd,s ) 
where pd,s is the disc and sphere densities, and Zd^s and 
Nd^s are the number of protons and the mean number 
of neutrons in the atoms of a disc and a sphere. The 
quantities pd,s are defined as pd,s = 'md,s/'mH where 
md,s are the mean masses of the disc and sphere atoms, 
and rriH is the mass of an atomic hydrogen, respectively. 
The values of Z/p and N/p for many elements with 
account of their isotopic composition can be found in 


The quantity I defined in ([3]) can be equivalently 
represented in the form [40] 

f °° ,/,,2 _ 1 

I{ma,D,Zi) = - du^ -e- 2 -aU.i ( 5 ) 

Jl U 


X (1 - . 


Using (ED and m, the magnitude of the differential ad¬ 
ditional force arising from the alternate interaction of 
a sphere with Au and Si sectors of the structured disc 
is equal to 


|‘AT’add(a)| = 


m^m 


-Cs{CAu - Csi) 


( 6 ) 


H 


du 


\/m 2 — 1 


( 1 - 


^ — 27nauD\ 


2 i?+a 

[i ?2 

Here, Cau and Csi are defined by ([1]) for Au and Si disc 
materials, respectively. 

It is convenient to introduce the new integration 
variable t = zi — a in the integral with respect to zi. 
After the integration with respect to t is performed, the 
differential force m is given by 

X [ (1 - d>{R, niau), 

Jl u 

where 

+ , 8 , 

Note that equation ED is the exact one. The differential 
force due to two-axion exchange can be obtained also 
using the proximity force approximation [52] , where the 
spherical surface is replaced with infinitesimally small 
plane plates parallel to the disc. This again results in 
ED, but with <!> = R. Such an approximate expression 
is only applicable under a condition R ^ 



In the next section, we apply the expression ED to 
calculate the differential force in the configuration of 
experiment [44] and obtain constraints on the coupling 
constants of axion-like particles to nucleons. 


3 Improved constraints on the coupling 
constants of axions to nucleons 


Now we take into account that in the experiment [33] 
the sphere was not homogeneous. It was made of sap¬ 
phire (AI 2 O 3 ) and covered with the layers of Cr and Au 
of thicknesses Z\cr = 10 nm and Z\au = 250 nm, respec¬ 
tively. Thus, the differential force between a sapphire 
core and a structured disc can be calculated by ED with 
Cs = CAI 2 O 3 where R is replaced with R — Aau — ^Cr- 
Then one should add to the obtained result the differ¬ 
ential forces between the structured plate and each of 
two spherical envelopes of external radia R — Aau and 
R of thicknesses Z\cr and Z\au, made of Cr and Au, re¬ 
spectively. These differential forces are calculated sim¬ 
ilarly using ED- For instance, the differential force due 
to a Au envelope is found by subtracting from ([7|) with 
Cs = Cau tfie differential force due to a Au sphere of 
radius R — Z\au placed at a separation a + Z\au from 
the structured plate: 

= 7. -V^^Au(Cau - Csi) (9) 


^ _ 1 
du A 


( 1 - 


— 2raauD\ 


[^{R, TUaU) - e maU)] . 

By adding up the differential forces from the sap¬ 
phire core and two spherical layers, we arrive at the 
following result valid in the experimental configuration 


133]: 


X [ X{maU), 

Jl u 

where 

X{z) = Cau [^R, z) - - Z\au, z)] (11) 

[d>iR-AAu,z) 

_g-2z4c,^(^_^Au-/\cr,^)] 


- 2 iAu - Acr,z). 

Note that the quantities C for Au, Cr, and AI 2 O 3 are 
defined by The values of Z/p and N/p for atoms of 

Au, Cr and Si, and for a molecule of AI 2 O 3 , as well as 
p for these materials, are presented in Table I of m- 
Now we are in a position to obtain constraints on 
the parameters of axion-like particles following from the 
measurement results of experiment [44] . Recall that nei¬ 
ther the Si substrate under the structured disc nor the 
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Fig. 1 Constraints on the coupling constants of axion-like 
particles to a proton or a neutron obtained from the Casimir- 
less experiment 1441 are shown as functions of the axion mass. 
The lines from bottom to top are plotted under the conditions 
alp = aln > alp, and glj, > gl^, respectively. The re¬ 
gions of the plane above each line are excluded and below 
each line are allowed. 

Au and Cr overlayers contribute to the measured dif¬ 
ferential force which is given by m and (HH). How¬ 
ever, when using the results of [23], it should be re¬ 
membered that the experimental separation distances 
between the rotating disc and the sphere are equal to 
z = a — dAu — d Cr = a — 160 nm. 

In the experiment [44] no differential force was ob¬ 
served. This means that the quantity m arising due 
to two-axion exchange between the test bodies was less 
than the minimum detectable force 

< SFia). (12) 

According to Fig. 3 of [44], at separation distances z = 
200, 400, 700, and 1000 nm (this corresponds to a = 
360, 560, 860, and 1160 nm) the minimum detectable 
force was equal to 6F — 0.2, 0.09, 0.12, and 0.12 fN, 
respectively. These values were determined at the 95% 
confidence level. We have found numerically the values 
of the axion-to-nucleon coupling constants gap, Qan and 
masses rria satisfying the inequality m with 
given by (nni) and (HU. The most strong constraints 
were obtained at a = 560 nm. 

The computational results for allowed and excluded 
values of the coupling constants 3ap(„)/(47r) as func¬ 
tions of the axion mass rua are presented in Fig. 1. The 
three lines from bottom to top are plotted under the 
conditions = gl^, gl^ » gl^, and gl^, » gl^, re¬ 
spectively. The regions of the plane {'ma,gap(n)) above 
each line are excluded by the experimental results, and 
the regions below each line are allowed. As can be seen 
in Fig. 1, with increasing TOq the strength of the ob¬ 
tained constraints quickly decreases, and they become 
not competitive. 


logic [5ap(n)/(47>-)] 



logic ["ia (eV)] 

Fig. 2 Constraints on the coupling constants of axion-like 
particles to a proton and a neutron obtained under the con¬ 
dition glp = gl^ from the magnetometer measurements 1341 
(line 1), from the Cavendish-type experiment I37II38I (line 2), 
in this work from the Casimir-less experiment 1441 (line 3), 
from measurements of the effective Casimir pressure I41II55L 
1561 (line 4), and from measurements of the lateral Casimir 
force between corrugated surfaces I42II57II58I (line 5). The re¬ 
gions of the plane above each line are excluded and below 
each line are allowed. 


It is interesting to compare the constraints of Fig. 1 
with previously obtained strongest laboratory constra¬ 
ints on axion-like particles. This comparison is pre¬ 
sented in Fig. 2 under the most reasonable condition 
gap = gan Eg. In Fig. 2, the line 1 shows the con¬ 
straints found from the magnetometer measurements 
using spin-polarized K and ^He atoms [34] (note also re¬ 
cent limits on a product of the pseudoscalar and scalar 
axion-to-nucleon interaction constants obtained |531l54j 
from the magnetometer experiment with ^He and ^^®Xe 
atoms). The line 2 shows the constraints derived [35] 
from the Cavendish-type experiment [37] . The line 3 
demonstrates the constraints obtained in this paper 
from the Casimir-less experiment m- This line repro¬ 
duces the lowest line of Fig. 1. The constraints found 
[41] from measurements of the effective Casimir pres¬ 
sure by means of a micromachined oscillator [551156] are 
presented by the line 4. Finally, the line 5 shows the con¬ 
straints derived [23] from measurements of the lateral 
Casimir force between sinusoidally corrugated surfaces 
[371(55] . The regions of (tos , 9ap(n )) plan® above all lines 
are excluded by the experimental results, and the re¬ 
gions below each line are allowed. 

As can be seen in Fig. 2, the obtained in this paper 
constraints of line 3 are significantly improved in com- 
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parison with the previously known ones in the wide 
region of axion masses from 1.7 x 10“^ eV to 0.9 eV. In 
the region from 1.7 x 10“^ eV to 4.9 x 10“^ eV our con¬ 
straints strengthen the constraints obtained |38] from 
the Cavendish-type experiment [87) whereas in the re¬ 
gion from 4.9 x 10“^ eV to 0.9 eV they are stronger than 
the constraints found m from measurements of the ef¬ 
fective Casimir pressure The largest strength¬ 

ening of previously known constraints by a factor of 60 
holds for TOq = 4.9 x 10“^ eV. Thus, the Casimir-less 
experiment leads to stronger constraints not only on 
the Yukawa-type corrections to Newton’s gravitational 
law [44], but on the coupling constants of axion-like 
particles to nucleons as well. 

4 Conclusions and discussion 

In the foregoing we have used the experimental results 
of recent Casimir-less experiment to derive the con¬ 
straints on the axion-to-nucleon coupling constants. The 
obtained constraints strengthen the previously known 
ones following from the Cavendish-type experiment and 
from measurements of the effective Casimir pressure. 
The strengthening up to a factor of 60 is achieved over 
the wide region of axion masses from 1.7 x 10“^ eV to 
0.9 eV. 

It should be remarked that the results of the Casimir- 
less experiment used here to obtain the constraints on 
axion-nucleon interaction are unambiguous in the sense 
that they do not depend on any theory of the Casimir 
force. This is different, for instance, from the constraints 
obtained |4T| from the measure of agreement between 
measured effective Casimir pressure and theory (see the 
line 4 in Fig. 2). It is well known that precise experi¬ 
ments on measuring the Casimir interaction between 
metallic surfaces agree with the extrapolation of a di¬ 
electric function to zero frequency using the plasma 
model, whereas a literally understood theory suggests 
to use the Drude model at low frequencies [5^1551156115^ 
[60ll61U62| . Recently the decisive experiment has been 
proposed [63| , where the theoretical predictions of both 
approaches differ by a factor of order 10^. The first 
measurements performed in the framework of this pro¬ 
posal are in favor of the plasma model extrapolation 
[M] . However, a fundamental understanding of phys¬ 
ical mechanisms behind this problem is still missing. 
Because of this, the obtained here constraints, which 
are independent on theory-experiment comparison for 
the Casimir forces, are of particular value. 

In the past, both measurements of the Casimir in¬ 
teraction and the Casimir-less experiment were used 
for obtaining constraints on the Yukawa-type correc¬ 
tions to Newtonian gravity [^I5^[551I551[551I56116711551 


155117(1] . In references [551155111711551155] and in this pa¬ 
per it is shown that the same experiments also lead 
to competitive constraints on the coupling constants of 
an axion-to-nucleon interaction. It should be remem¬ 
bered, however, that the Yukawa potential arises due 
to exchange of one scalar particle, whereas the spin- 
independent potential dm, used to constrain the param¬ 
eters of axion-like particles, results from the exchange of 
two particles. The latter makes the obtained constraints 
relatively weaker. Thus, in the future it seems promis¬ 
ing to perform measurements of the Casimir interac¬ 
tion and Casimir-less experiment using the polarized 
(magnetized) test bodies. In this way one could obtain 
much stronger constraints on the parameters of axion- 
like particles by exploiting spin-dependent interaction 
potential arising due to one-particle exchange between 
nucleons. 

Acknowledgements The authors are grateful to R. S. Decca 
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